Discontinuous current-phase relations in small ID Josephson junction arrays 
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We study the Josephson effect in small one-dimensional (ID) Josephson junction arrays. For weak 
Josephson tunneling, topologically different regions in the charge-stability diagram generate distinct 
current-phase relationships (I$s). We present results for a three-junction system in the vicinity of 
charge degeneracy lines and triple points. We explain the generalization to larger arrays, show that 
discontinuities of the 1$ at phase n persist and that, at maximum degeneracy, the problem can be 
mapped to a tight-binding model providing analytical results for arbitrary system size. 

PACS numbers: 73.23.-b, 74.50.+r, 74.81.Fa 



Introduction. — The Josephson current-phase relation 
(I<f>) contains information about the microscopic nature 
of the Cooper pair transfer between superconductors. 
Specifically, systems more complicated than a single tun- 
nel junction generally exhibit non-sinusoidal I<J>s [l[. A 
pertinent example consists of the double junction formed 
by one superconducting grain coupled to two leads 0, Q . 
In this Letter, we go beyond the single-grain case and 
investigate the Josephson effect in small ID Josephson 
junction arrays (JJAs), cf. Fig. [Ha). 

As demonstrated for the two-junction case, charging 
effects can play an important role in the transport of 
Cooper pairs due to the competition between Coulomb 
blockade and tunneling of charge. This can significantly 
modify the I<& and the critical current 0, To be 
precise, the 1$ of the two-junction system differs max- 
imally from a sinusoid at the charge degeneracy point. 
Given that both junctions are equally strong, the I<I> is in 
fact discontinuous at <fi = it, where cj> denotes the gauge- 
invariant phase difference across the entire array. 

These results are obtained in the particularly inter- 
esting regime where superconductivity is fully developed 
(T <C Ec <C A, the respective quantities denoting the 
temperature, relevant charging energy, and supercon- 
ducting gap), but charging effects are strong compared to 
Josephson tunneling, Ec 3> Ej. In this case, the chargc- 
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FIG. 1: (Color online) (a) Progression of small ID Joseph- 
son junction arrays with increasing number of junctions N. 
(b) Detailed schematic of the three-junction JJA. Electronic 
transport across the junctions is due to inter-dot (H m ) and 
lead-dot (Hl,r) tunneling. The electrochemical potentials of 
the dots can be adjusted by separate gate electrodes. 



stability diagram [3] is an appropriate starting point for 
identifying the relevant states, and for classifying the dif- 
ferent types of I$s. For instance, the charge-stability 
diagram of the two-junction system consists of isolated 
points along the gate- voltage axis, specifiying the poten- 
tials where two charge states become energetically degen- 
erate. In the vicinity of these charge-degeneracy points, 
two charge states participate in the Cooper pair transfer 
and the I<I> is strongly non-sinusoidal. Away from the 
degeneracy points, the system remains in a single charge 
state and exhibits a sinusoidal 1$. 

This reasoning can be extended to ID JJAs with larger 
numbers of grains, making the connection between the 
limiting cases of a single grain and the limit of an in- 
finite array. It is important to note, however, that the 
regime considered here is different from the self-charging 
and nearest- neighbor models for infinite ID JJAs [H, 0] • 
In our case, the total number of junctions in the array 
determines the maximum number of charge states which 
can become degenerate. The maximal degeneracy is gen- 
erally represented by a point in the charge-stability dia- 
gram. Lower-order degeneracies are associated with hy- 
pcrplanes of different dimensionalities (lines and planes 
being simple examples), generating distinct I$s. The ac- 
cessibility of such degeneracies up to fourth order has re- 
cently been demonstrated in normal-state triple dots 0] . 
We show that discontinuities in 1$ persist beyond the 
single-grain case. We present a detailed analysis of the 
three-junction system shown in Fig. [ljb), and then de- 
scribe the generalization to larger systems; at maximum 
degeneracy we obtain a useful tight-binding approach. 

Model. — For three junctions, the system consists of 
two superconducting quantum dots a and b in series 
with two macroscopic superconducting leads L and R, 
closely resembling devices previously studied in various 
contexts in experiments @, (3 EH- Th e dots and leads 
are coupled through Josephson junctions with Joseph- 
son energies Ej a and junction capacitances C a , where 
a = L,R,m signifies the left, right, and middle junction. 
Additionally, the two dots are capacitively coupled to 
separate gate electrodes, enabling the individual tuning 
of their electrochemical potentials. 
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FIG. 2: (color online) (a) Lumped element circuit for the su- 
perconducting double quantum dot. The supercurrent in the 
loop is induced by applying an external magnetic flux $ and 
can be measured by an inductively coupled tank circuit (not 
shown) . (b) Charge-stability diagram for the superconducting 
double-dot system. The dashed regions mark the regimes (i)- 
(iii) away from all degeneracies, close to a charge-degeneracy 
line, and close to a triple point of the honeycomb pattern. 



We focus on the scenario where the superconductiv- 
ity is fully developed, captured by BCS theory, and not 
impeded by charging effects (T <C Eq <C A). In this 
case, breaking of Cooper pairs due to thermal fluctu- 
ations or charging effects is negligible, and the system 
Hamiltonian can be derived in a convenient way via cir- 
cuit quantization cf. Fig. [UJa). The Hamiltonian 
consists of charging and Josephson tunneling contribu- 
tions, H = H c h + ^Z a H a . The first term describes the 
charging energy of the double-dot system 13, , 



H ch = E Cj(nj - n gjf + E m 17 (hj - n gj ), (1) 

j—a,b j—a^b 

where h a , fit, count the number of extra Cooper pairs 
on dots a and 6, respectively. For convenience, we have 
introduced the abbreviations 
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for the partial charging energies. The tunable offset 
charges are denoted by n ga = C ga V ga /2e, and C^ a ,b = 
C m + Cl.r + C ga ,b- In the following, we will restrict to 
the case where all junction capacitances are finite. The 
corresponding charge-stability diagram is shown in Fig. 
IHh). For finite C™, it forms a honeycomb pattern typi- 
cal of double dots [4(], and features three distinct regimes: 
(i) regions away from any charge degeneracies, (ii) charge 
degeneracy lines, and (iii) triple points. 

The tunneling of Cooper pairs is captured by the 
Josephson terms H a , where 



H L = -E JL cos(ip L - tp a ), H R 
H m = -E Jm cos(ip a - tp b ). 



-E JR C0S(V3R - ip b ), 
(3) 



Here, the quantities tpj with j = a, b, L, R denote the 
phases used to build the gauge- invariant phase differences 



across the junctions. Note that the leads L and R are as- 
sumed to be macroscopic superconductors with negligible 
charging energy. As a result, as opposed to the phases 
on the grains their phases may be treated "classically" . 

Current-phase relations. — The I*!* may be probed in a 
loop configuration as shown in Fig. [DJa), employing the 
"rf technique" discussed in Ref. Flux quantization 
relates the external flux <E> to the phase difference <f> via 
4> — <Pl — <Pr — 27r$/$o, where $o = h/2e denotes the 
flux quantum. Based on the above model, we evaluate 
the I$s in the regimes (i)-(iii) determined by the charge- 
stability diagram. For (i), i.e. away from charge degen- 
eracies, it is straightforward to see that each Cooper pair 
must coherently tunnel from one lead to the other, occu- 
pying the grains only in virtual intermediate states. The 
result is a sinusoidal I4>, analogous to the single-grain 
case away from charge degeneracy. Here, however, the 
critical current is suppressed by a factor (Ej a /Ec) 2 - 

Close to charge-degeneracy lines [regime (ii)], two 
charge states participate in the Josephson tunneling. 
Without loss of generality, we may focus on the gate- 
voltage region marked in Fig. [21b), where the charge 
states (1,0) and (0,1) with an additional Cooper pair 
on either dot a or 6 are close to degeneracy. Our 
assumption of weak Josephson tunneling then requires 
Ej a < Ec = iam n {E„ - E i,E n - E w }. Here, E jk 
stands for the charging energy of the state with j (k) ad- 
ditional Cooper pairs on dot a (b), and n enumerates all 
charges states different from (0,1) and (1,0). The cor- 
responding effective Hamiltonian, acting in the reduced 
Hilbert space spanned by the charge states (0, 1) and 
(1,0), can be constructed systematically order for or- 
der, cf. Carrying out this procedure up to order 
0(Ej a / Ec) and discarding of phase-independent terms 
merely renormalizing the charging energies, we obtain 

Here, we have employed a pseudospin description for the 
charge states (1,0) and (0, 1), reinterpreting them as the 
spin-up and spin-down eigenstates of the operator a z . 
Any deviation from charge degeneracy (tunable by the 
gate voltages) leads to a "magnetic field" ho = Eiq—Eqi. 

Transfer of Cooper pairs between dots and leads is de- 
scribed by the terms ~ J, J' in Eq. Q. Here, J' origi- 
nates from inter-dot tunneling processes, Fig. [3)^a), and 
J is generated by dot-lead tunneling, Fig. [31(b). In the 
spirit of the nomenclature for single electron tunneling 
processes, we may refer to these different processes as se- 
quential tunneling and cotunneling of Cooper pairs. The 
energies J, J' > are given by 
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Ejj - E w 



(5) 



and J' = Ej m /2. Note that the ratio J/J' weakly de- 
pends on the gate voltages. Additional tunability of this 
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FIG. 3: (color online) (a,b) Leading-order processes relevant 
to charge transfer in the three-junction system, close to a 
charge-degeneracy line. Depicted are the processes starting 
in the initial charge state (0, 1). The inter-dot tunneling (a) 
contributes to J', the dot-lead processes (b) contribute to J. 
(c) Energies E±(cj>), Eq. (|6j), of the two-level system close to 
a charge-degeneracy line. The eigenenergies are plotted as 
functions of phase for the resonant case J = J'. The solid 
(dashed) lines correspond to ho = (J/2), the arrow marks 
the position of the low-temperature 1$ discontinuity at n. 



parameter may be achieved by terminating the super- 
conducting leads in a 2D electron gas, cf. e.g. 



151 ] . which 



enables the direct tuning of tunneling strengths. 

By diagonalization of Eq. ^ we obtain the eigenener- 
gies plotted in Fig.^c), 



E ± (c/>) = ±y/(h /2)*+e 2 JJI (<l>), 



(6) 



where £jj'(</>) = [J 2 + (J 1 ) 2 + 2JJ' cos0] 1/2 . Note that 
for ho = the two states become degenerate at <f> = tt. 
In the Heisenbcrg picture, the current is 



2ie 

-h J{a 1 



(7) 



With this, the supercurrent can be computed for both 
zero and finite temperatures. For T = 0, it is given by the 
ground state expectation value of the current operator, 



2c 



J J 1 sin ( 



(8) 



Alternatively, the zero-temperature supercurrent may be 
obtained via Io{4>) = xi^T' lading to the same result. 
At finite temperatures small compared to A, quasiparti- 
cle excitations may be neglected and the 1$ is obtained 
by averaging with respect to the equilibrium density ma- 
trix, resulting in I(4>) = I ((f>) tanh \J3E-(<j>)]. This 1$ is 
depicted for T = in Fig. Ufa). Note that due to the gen- 
eral symmetry properties of I<E»s [l[ it is sufficient to plot 
the range < < 7r. This 1$ is intimately related to the 
result found in the single-grain case at charge-degeneracy 
points 0]. In both cases, Cooper pair transfer involves 
the switching between two charge states - occupation and 
de-occupation in the single-grain case, occupation of the 
right or the left grain in the present case. This analogy 
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FIG. 4: (color online) Current-phase relations for the three- 
junction JJA at zero temperature, (a) close to a charge- 
degeneracy line, (b) at a triple point, (a) 1$ for different 
ratios J' /J of the inter-dot and dot-lead coupling, (b) 1$ for 
Ej = Ejl = Ejr and different ratios Ej m /Ej. 



and the underlying pseudospin description are the ori- 
gin for the I^s to be structurally identical. In particular, 
in close resemblance to the discontinuity for the single- 
grain system, we find that, given J = J', the 1$ of the 
double-grain system undergoes a discontinuity at <fi = tt, 



I (4>) = ^Jsm(4>/2) for 
n 



TT < < TT, 



(9) 



which gets broadened by temperature. However, simi- 
lar to the situation of the sinusoidal relation in regime 
(i), the critical current at the charge-degeneracy line re- 
mains suppressed by a factor Ej a /Ec, see Eq. (|SJ). It is 
interesting to note that similar 1$ discontinuities are also 
present in classical and ballistic point contacts [l], [lj], 17 1 . 



A few comments on the 1$ © and its crossover to 
a sinusoid arc in order. First, we remark that Eq. ([8]) 
is invariant under interchange J J', implying a one- 
to-one correspondence between the supercurrent through 
a double-dot system with large inter-dot tunneling but 
small dot-lead coupling, and the current for small inter- 
dot tunneling but large dot-lead coupling. Second, the 
effect of detuning h from the charge-degeneracy line can 
always be mimicked by adjusting the parameters J and 
J', so that the consideration of the ho = case is suffi- 
cient. Third, for J'/J^l, we asymptotically recover a 
sinusoidal 1$, 

2e 

I(4>)~— min{ J, J'} sin tanh \(3 max{ J, J'}] . (10) 

h 1 J 

This is in accordance with the behavior of the single-grain 
1$ when making the Josephson energies in the left and 
right junction very different 

In the regime (iii) , the system is tuned to a triple point 
of the charge-stability diagram, where three charge states 
become degenerate. We focus on the triple point spec- 
ified by a circle in Fig. [2Kb) , where the corresponding 
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Hamiltonian may be represented as 



1 



2E, 



DO 



-E. 1L e^ -E. 
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2E i 



(11) 



For simplicity, we present an analytical 1$ expression 
only for T = 0, exact degeneracy, and identical junc- 
tion strengths. Denoting Ej = Ejl = Ejr = Ej m , one 
finds 



2e E, 



-sin((/>/3) for — tt < <p < tt, 



(12) 



giving the maximum possible supercurrcnt in the three- 
junction device. Finite temperatures, detuning from the 
triple point, and differing junction strengths generally 
lead to a broadened 1$ with more complicated expres- 
sions. Eq. (fT2"|) . along with examples for the situation of 
differing junction strengths, is plotted in Fig.[3Jh). 

Generalizations. — For an array consisting of N junc- 
tions, maximally N states can become degenerate at iso- 
lated points in the charge stability diagrams. The I<I> 
hierarchy then consists of N different I<E»s whose critical 
currents differ in powers of the small parameter E j / Ec ■ 
Away from degeneracies, we find that the 1$ is sinusoidal 
with critical current suppressed by a factor (Ej / Ec) N ~ x ■ 
In the vicinity of two-fold degeneracies, the 1$ takes the 
form Eq. ([9|) with discontinuity at phase tt; the criti- 
cal current in this region is suppressed by (Ej / Ec) N ~ 2 ■ 
When increasing the degree of degeneracy, this hierarchy 
continues and discontinuities persist. 

On the uppermost hierarchy level, the critical current 
is not suppressed and the discontinuity is fully developed 
when all junctions have the same Josephson energy Ej. 
In this case, the system can be mapped to an effective 
tight-binding model with periodic boundary conditions, 



tb 



4e 



e -i<t>/N I v . 



i)M 



H.c. 



(13) 



Here, the states | v ) denote the N degenerate charge 
states and can be pictured as the states with one extra 
Cooper pair on one of the grains or the lead. The eigenen- 
ergies of the system are given by e n = —Ej cos(k n + (f)/N) 
with k n = 2irn/N and n = 0, • • ■ , (N — 1). The resulting 
zero-temperature 1$ is 



r (Y> 



(0) = y^sim>/7V) for -tt< 



< 7T, 



(14) 



which generalizes Eq. (|12[) for arbitrary number of junc- 
tions TV, and approaches a shallow sawtooth for large N. 
Owing to the general 2tt periodicity of I<E>s [l[ , Eq. (fH)) 
proves that discontinuities of the 1$ at tt persist. It is 
interesting to note that Eq. lfl"4|) is reminiscent of the 
classical limit for a ID J J A 0]. The crucial difference is 
that, in the present case, there is a probability 1/7V to 



occupy any of the charge states. As a result, the critical 
current scales with N~ 2 instead of ./V -1 for large N. 

Conclusions. — We have studied I<f>s in finite-size ID 
JJAs. In the limit Eq 3> Ej, the charge-stability di- 
agram serves as an important tool to identify different 
types of I<E>s. We have established a hierarchy of I<I>s, 
where each hierarchy level belongs to a different degree 
of charge degeneracy and has a distinct magnitude of the 
critical current. Under specific conditions and zero tem- 
perature, these I<f>s display discontinuities at phase tt, 
originating from a crossing of the lowest two eigenener- 
gies. Maximum charge degeneracy allows for a mapping 
to a tight-binding model, and results in a peculiar I<I> 
converging to a sawtooth function with a critical current 
scaling with N~ 2 . With experiments on small [H, QlJ 111 
and large arrays [l8| having been demonstrated, we be- 
lieve that our predictions will have an experimental im- 
pact in the near future. 
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